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Abstract 

Recent data from heavy ion collisions at RHIC show unexpectedly large near-angle correlations that broaden 
longitudinally with centrality. The amplitude of this ridge-like correlation rises rapidly with centrality, reaches a max- 
imum, and then falls in the most central collisions. In this talk we explain how this behavior can be easily understood 
in a picture where final momentum-space correlations are driven by initial coordinate space density fluctuations. We 
propose vl/el part as a useful way to study these effects and explain what it tells us about the collision dynamics. 
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Introduction: The motivation for the construction of the Relativistic Heavy Ion Collider at Brookhaven National 
Laboratory was to collide heavy nuclei in order to form a state of matter called the Quark Gluon Plasma (QGP) HI . 
These collisions deposit many TeV into a region the size of a nucleus. The matter left behind in that region is so hot 
and dense that hadronic matter undergoes a phase transition into a form of matter where quarks and gluons are the 
relevant degrees of freedom, not hadrons |2|. This is the state of matter that existed when the universe was less than a 
microsecond old and still very hot. 

Correlations and fluctuations are an invaluable tool for probing the dynamics of heavy-ion collisions. Data from 
the experiments at RHIC reveal interesting features in the two-particle correlation landscape J3]|4). Specifically, it 
has been found that correlation structures exist that are unique to Nucleus-Nucleus collisions. While two-particle 
correlations in p+p and d+Au collisions show a peak narrow in azimuth and rapidity, the near-side peak in Au+Au 
collisions broadens substantially in the longitudinal direction and narrows in azimuth. An analysis of the width of the 
peak for particles of all pj finds the correlation extends across nearly 2 units of pseudo-rapidity (A77 = 2) (4). When 
triggering on higher momentum particles (pj > 2 GeV/c for example), the correlation extends beyond the acceptance 
of the STAR detector (A/7 < 2) and perhaps as far as A77 = 4 as indicated by PHOBOS data |4). 

STAR has found that the amplitude of this correlation shows a rather rapid rise with collision centrality [3 1 before 
reaching a maximum and falling off in the most central bins. This drop in the most central bins shows up for both 
■Js m = 200 and 62.4 GeV but is often overlooked. In this talk we present a geometric explanation for the centrality 
dependence of the ridge amplitude. We'll use the centrality dependence of V2/S2, dN/dy, and the third harmonic 
participant eccentricity e^ part to predict the amplitude of the near-side ridge correlation (Aj). Given the apparent 
relevance of initial-state density fluctuations to the final, momentum-space correlations, we advocate the transfer 
function vl/s^ pan as a valuable observable for studying the length scales in heavy-ion collisions. We extract the 
transfer functions from intermediate pj di-hadron correlation data where evidence has been presented for conical 
emission Q- 

Three Premises: It has been shown by the STAR collaboration that the second harmonic component of the near- 
side ridge in the two particle correlations can account for nearly all of the difference between the two and four particle 
cumulant V2 Q. It was argued that vi fluctuations must therefore be tiny: many times smaller than the fluctuations 
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predicted from eccentricity models (S)- In that case, a major revision of our understanding of heavy-ion collisions 
would be required. We've argued previously, however, that eccentricity fluctuations can give rise to v„ fluctuations for 
more harmonics than just n — 2, and that those fluctuations could therefore be the source of the near-side ridge J9] 
(especially if the v„ fluctuations depend on pseudo-rapidity difference A77 i.e. (v n (r])v„(ri + Arf)) = /(A77)). This idea 
has been born out by calculations from several groups IfTUl 151 [TTl [T2l . To test this conjecture, we will attempt to 
explain the centrality dependence of the near-side ridge amplitude A\ from eccentricity fluctuations, we start with 
three simple premises: 

• the expansion of the fireball created in heavy-ion collisions converts anisotropies from coordinate-space into 
momentum-space, 

• the conversion efficiency increases with density, 

• and the relevant expansion plane is the participant plane. 

The participant plane can be defined for any harmonic number and a system with a lumpy initial energy density 
will give rise to finite participant eccentricity at several harmonics |5][T3]. This becomes conceptually clear when 
eccentricity is recast in terms of a harmonic decomposition of the azimuthal dependence of the initial density. Our 
calculation of the centrality dependence of A\ will depend on the higher harmonic terms in an eccentricity model. 

Higher Harmonics, Even the Odd: Fig. Ill (left) shows the n'^-harmonic participant eccentricity {sf u , art ) as de- 
fined in Ref. [5 1 for central Au+Au collisions from a Monte-Carlo Glauber model. Typically the participant eccentric- 
ity is calculated based on the positions of point-like participants (that is the participant is said to exist at a precise x and 
y position). Those results are labeled r par , - 0.0 fm. One can also calculate the eccentricity from a more physically 
realistic model with participants smeared over some region. This is done by treating each participant as many points 
distributed within a disk of some finite radius r part . Increasing r part washes out the higher {s^ part ) terms. The right 
panel of Fig. [T] shows the ratio of {s^ part ) for a given r pan value divided by {sf t parl ) for r part = 0. The curves are 
labeled l m f p instead of r pan and the ratio l m f p /idea\ for reasons explained below. 




n n 

Figure 1: Left panel: n ) for 0-5% central Au+Au collisions from a Glauber Monte Carlo where participants are either treated as point-like or 
they are smeared over a region of size r pml . We include the oblate shape of Au nuclei in our calculation. Right panel: the ratio of the curves on the 
left to the point-like curve. 

In this calculation we introduced the length scale r parl causing the higher terms in (e% part ) to be washed out. That 
effect is more general though and we believe it is important for understanding correlations and v„ fluctuations. One 
can also consider what happens when particles free-stream for some amount of time T/ s before they interact; that 
will also introduce a length scale ctj s which leads to a similar reduction of higher terms IfTZl . One can also consider 
the effect of a mean-free-path (l m f p ) on the ability of the fireball to convert higher {e\ part ) terms into v\ |[T4l[T5l . If 
a particle on average travels for a distance l m f p between interactions, it is clear that higher (sj; part ) terms will not 
become manifest in v\. Fig. [2] shows a schematic illustration of this idea. If our probe has a l m / p in the fireball, then 
we will be blind to features smaller than l m f p . All these effects will act to wash out the higher harmonics so we expect 
v\ to drop with n. Since v\ is related to dN/dA(p by a Fourier transform, and since a Fourier transform of a Gaussian 
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is a Gaussian, all we need to reproduce the near-side Gaussian ridge in two particle correlations, is for v\ to drop with 
n with an approximately Gaussian shape. The right panel of Fig. [T] shows that this is a reasonable expectation from 
viscous effects. 



Higher harmonics probes smaller length-scales. 




n=2 n=3 n=4 n= 1 n=IS 



Figure 2: Schematic illustration of the interplay between a length scale like the mean-free-path and the higher harmonics. 

The ridge amplitude A\ is found by measuring hp/ ^Jp re f (the pair density p minus the reference pair density p re f 
scaled by ^p re f) vs h<p and A77 ifTBI . A fit function is devised to describe the correlation. The fit function has a h<p 
independent term, a cos(A^) and cos(2A0) term, and a near-side 2-D Gaussian. The fit function describes the data 
well 1161 . Here we work with the conjecture that the 2-D Gaussian is a manifestation of (s^ part ). Based on this 
conjecture, we can try to calculate the centrality dependence of A 1 . Our result for A 1 will be related to v\ so we need 
to know the conversion efficiency c of (sl pan ) into v\. The conversion efficiency will depend on particle density. In 
Fig. [3] we show v^/e vs density (l/S)dN/dy from the STAR Collaboration |17| with a fit function to parameterize 
the data. The figure shows V2{4) divided by the eccentricity calculated with respect to the reaction plane from a CGC 
model Ifl8l . We will take this to estimate our conversion efficiency with the understanding that the different results 
obtained from CGC and Glauber models gives rise to at least a 30% uncertainty in the correct conversion efficiency. 




Figure 3: The ratio of viW over £ stl f from a CGC model vs {} I S )dN j dy. 



Since a cos( h(/>) and cos(2A0) term have already been subtracted from hp/ -\Jp re f in order to obtain A 1 , our estimate 
of A\ can be made simpler if we predict the n — 3 component of the near-side ridge and scale that up to get the full 
amplitude. The n — 3 component is found from 

— f -^-cos(3A0)rfA0 = O.O39Ai. (1) 

The azimuthal width <x of the near-side Gaussian is weakly dependent on centrality so we use a typical value of 
cr = 0.65. The factor of 0.039 will be used to relate our prediction for (s^ part ) to the amplitude A\, Since hp/ yjPref 
is a per-particle measure instead of a per-pair measure, we need to include the particle density po = j~ ■ Putting 



/ Nuclear Physics A 00 (2011) 1^6\ 



4 



(gj ) together with the conversion efficiency c, particle density po, and the factor of 0.039 to map from the n — 3 
component to a Gaussian amplitude, we find 

Ai*Poc4 part /0.039. (2) 

We take po from data, c from Fig. pi and (s^ ) from our Monte-Carlo Glauber model. 

The Rise and Fall of the Ridge: The right panel of Fig. [4] shows our estimate of the ridge amplitude A\ based 
on s^ part vs centrality parameter v = 2Nt,i„/N parr . The left panel shows the preliminary STAR data. Our estimate for 
the amplitude is only approximate since we've applied the conversion efficiency for n — 2 to n — 3 and as discussed 
before, we expect the efficiency to drop with n. This will lead us to overestimate the contribution from {e 2 part ) to 
the ridge. On the other hand, we used a CGC based eccentricity model to extract the conversion efficiency but a 
Glauber model for {e\ parl ). This should cause us to under-estimate the ridge contribution; to some extent canceling 
the previous overestimate. We find that our estimate of A\ agrees quite well with data. The centrality dependence is 
particularly interesting: our A\, like the data, starts at a small value and rises much faster than expectations from a 
Glauber Linear Superposition model (GLM) which assumes that correlations grow as Nbi n /N part . The rise continues 
until A\ reaches a maximum near v = 5, then A\ falls again. This rise and fall is also seen in the preliminary 62.4 
GeV data [16| and has no natural explanation in any of the other proposed scenarios for the ridge formation. In our 
picture, the rise and fall is related simply to the geometry and it's fluctuations. {e\ pan ) falls with N part since the larger 
N part leads to smaller fluctuations. But N par ,(s^ ) first rises then falls. This rise and fall is due the asymmetry of 
the overlap region which shows up even for (s^ ). Since both c and po are increasing with centrality, the product 
of pqc{e\ pan ) rises until very central collisions and then falls as shown in the figure. The observation that the rise and 
fall shows up in the near-side ridge amplitude suggests that the near-side ridge is likely dominated by initial geometry 
fluctuations. 

V 0_8r 1 




Figure 4: Left panel: measured near-side ridge amplitude A[ vs centrality measure v = 2Ni,i,JN part . Right panel: contribution to the ridge from 
initial state geometry fluctuations. 

Transfer Functions: Our explanation for the centrality dependence of A; is based on the three simple premises 
listed earlier and provides a natural explanation for the rise and fall of the ridge. Initial estimates of the amplitude agree 
to within our uncertainties. This suggests that geometry fluctuations in the initial overlap region are converted into 
momentum space giving rise to the near-side ridge structure. We discussed that we expect the conversion efficiency to 
drop with n since effects like initial-state free-streaming and mean-free-path will wash out the higher harmonic terms. 
Measuring the conversion efficiency c„ = - ''"* 2 | 2 p as a function of n, centrality, and particle kinematics will provide 
information on those effects. It will be particularly interesting to measure c„ as a function of Arj to understand how 
de-coherence affects manifest in the longitudinal direction. 

In Fig.plwe show the conversion efficiency v„{2} 2 / part for intermediate pj di-hadron correlations from STAR J6|. 
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Figure 5: The transfer function v„{2| 2 /£^ p(]rf for intermediate pr di-hadron correlations from 20-60% central 200 GeV Au+Au collisions |61 . 

The panels show various combinations of trigger particle and associated particle pj selections. Intermediate pj di- 
hadron data rather than showing an away-side Gaussian, show two peaks shifted to either side of A(f> = n. It has 
been suggested that this is evidence of conical emission on the away-side of the higher pj trigger particle Ifl9l l20l . 
The previously ignored effects of v? which were, however, could explain some or all of these novel structures |9|. 
To investigate these structures, we look at the pj dependence of c n (p'j S , p a T ssoc ). By plotting v„{2} 2 / si pan from these 
correlations, we hope to see what portion of those correlations can be explained by geometric effects. For relatively 
lower momentum cuts, we find that c„(2.5 GeV, 0.5 GeV) behaves much as we expect from mean-free path effects 
for example; the higher terms drop off monotonically similar to Fig. [T| (right). But for higher p a T ssoc cuts c n shows a 
pronounced peak at n = 3. See c„(2.5 GeV, 2.5 GeV) for example. This feature disappears again however when a 
large enough pj g cut is applied. For p'j S > 6 GeV for example, we see only effects due to the presence of correlations 
from jets. It remains interesting to speculate about the possible source of that local maximum at n = 3 for interme- 
diate p T di-hadron correlations (whether it's due to a suppression of the lower harmonic super-horizon modes dl . 
mach-cones llZOI or some other acoustic effects [22]). It seems to be larger than trivially expected from el part . We 
note however that a complete investigation of the systematic errors on our Glauber Model has not been carried out. 
We include the oblate shape measured for the Au nucleus. That oblateness leads to a large enhancement of the n — 2 
component of £ npan . Reducing the oblateness of the Au nucleus could therefore reduce or eliminate the prominence 
of the n — 3 term in c„. This remains for further investigation. 

Conclusions: We presented the participant eccentricity vs harmonic when the participants are treated as point-like 
or smeared over a radius r pan . The larger values of r part wash out the higher harmonic eccentricities. We argued that, 
similarly, a large mean-free-path should wash out higher harmonics of v„. Such an effect could lead to a Gaussian 
peak in two particle correlations at A<f> = 0. We've calculated the contribution to the near-side Gaussian peak that we 
expect from initial density fluctuations. We based our calculation on three premises 1) that the expansion of the fireball 
created in heavy-ion collisions converts anisotropics from coordinate-space into momentum-space, 2) the conversion 
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efficiency depends on particle density, and 3) the relevant expansion plane is the participant plane. Following these 
premises, we find that the near-side peak from density fluctuations should rise rapidly, reach a maximum just before 
the most central events, then fall. Our estimate of the magnitude is in agreement within our uncertainties with the 
available data and the shape matches that seen in the data. This is the only calculation to correctly describe the rise and 
fall of the ridge amplitude. We conclude therefore that density fluctuations are likely the dominant source for the low 
Pt ridge-like correlations. We have also shown the ratio of the final momentum space correlations v„{2} 2 to the initial 
coordinate-space eccentricities e\ part for intermediate pj di-hadrons. We find that for intermediate pj correlations 
the n — 3 term is larger than n — 2. This suggests that even after taking into account initial density fluctuations, some 
interesting signal may exist in the intermediate p T di-hadron data. 
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